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Abstract 

The properties of the evolution equation have been analyzed. The 
uniqueness and the existence of solution for the evolution equation 
with special value of parameter characterizing intensity of change of 
external conditions, of the corresponding iterated equation have been 
established. On the base of these facts taking into account some prop- 
erties of behavior of solution the uniqueness of the equation appeared 
in the theory of homogeneous nucleation has been established. The 
equivalence of auxiliary problem and the real problem is shown. 



1 Formulation of the problem 

In P the evolution equation for the nucleation period has been derived. Here 
we shall analyze the properties of the evolution equation and the uniqueness 
and existence of the solution of special problems appeared in the theory of 
nucleation. 

Consider the initial problem formulated in pQ 
• The following equation is given 

g{z) — f* ( ( z — x) 3 exp(cx — Tg(x)/&*)dx (1) 

J — oo 

with positive parameters T, c, which are chosen to 
satisfy condition 

^l~ = <*./r (2) 
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Equation (|TJ) and condition (J2J) form the initial problem. 
This problem as well as all other ones is considered in a class of continuous 
functions. 

Instead of the function g we introduce the function G = Tg/Q*. Then we 

get 

G(z) = (/„r/$„) f (z - x) 3 exp(cx - G(x))dx 

J —oo 

dG(z) . 

—r \Z = = C 

dz 

We change the variables 

z - (/,r/$,) -1/ V 
x - (/ # r/^)-vv 

The new variables will be marked by the same letters z, x. 
Then we get 

G(z)= f* {z-x) 3 exp{cU'S/^*y 1/i x-G{x))dx 

^U=o = c(/.r/$.)" 1/4 

az 

We denote c(/ ;(! r/$ ;)! )~ 1 / 4 also as c. The value G will be denoted as g. We 
come to the following problem 

• The following equation 

g(z) — / (z — x) 3 exp(cx — g{x))dx (3) 

J — oo 

is given with the positive parameter c, chosen to satisfy 

dg{z) 



dz 



(4) 



The equation Q and the condition (jJJ) form the reduced problem. 

So, here only one parameter c remains. 

It is more convenient to consider instead of g the function 

4> = cx — g 

Then we come to the following problem (Problem A) 
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• The following equation 



cz = 



4>{z) + I (z — x) 3 exp((p(x))dx 



(5) 



— oo 



is given with the positive parameter c, chosen to satisfy 



Equation (j3j) and the condition © form the problem A. This problem 
will be the auxiliary problem. 
Consider now the problem B 

• The following equation 



is given 

Only the equation Q forms the problem B. 
We see that the problem B has no parameters. 

The main goal of investigation will be to see the uniqueness of the solution 
of the problem B. 

2 Solution of equation (G3) with fixed c 

2.1 Some properties of solution 

At first we consider (J3J) with some positive c. 
We rewrite (0) for function g in the form 



3 / y 2 exp((j)(y))dyz = <p{z) + / (z — x) 3 exp(<fi(x))dx 





(7) 




One can see the following property 
If the solution exists, then 



g{z) > 



for any z 
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This property goes from the positivity of sub-integral function. 
We introduce the iterated equation 

/Z rx 
(z - x) 3 exp(cx - (x- yf exp(cy - g(y))dy)dx (8) 
-oo J — oo 

We also introduce the nonlinear operator P according to 
P( s ) = ( z ~~ x ) 2 exp(ca; — s(x))dx 

J —oo 

One can prove the following fact: 

If for any x the following inequality 

9i(x) > g 2 {x) 

takes place, then 

P( 9l (x)) < P(g 2 (x)) 

takes place for any x. 

It follows from the explicit form of P. 
One can also see the following fact 
If the solution exists then 

g(z) < exp(c^)/6 

for any z 

It follows from the positivity of exponent in the sub-integral function. 
We construct iterations according to 

9i+i = P{9i) 
As the initial approximation we choose 

2o = 

We see that 

92 > 

for any value of the argument 

It follows from the positivity of exponent and the whole sub-integral func- 
tion. 
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Form the statements formulated above it follows the following chain of 
inequalities 

= go < 92 < 9a < •••• < 9 < •••• < 95 < 93 < 9i 

for any value of the argument. 

So, the odd iterations converge and the even iterations converge. So, the 
solution of the iterated equation exists. But the uniqueness isn't yet proven. 

For any initial approximation which is always positive the iterations with 
initial go = estimate them from above. 

Since the r.h.s. is positive and solution has to be positive, the iterations 
with any initial approximation have to be positive starting from some number 
(the asymptotics at small arguments are written explicitly) . Correspondingly, 
later the iterations will be estimated by the already constructed iterations. 
So, the convergence can be estimated to convergence of constructed itera- 
tions. 

Remark: The analogous iterations have been constructed by F.M.Kuni 
in 1984, but for the problem B. For the problem B the monotonic functions 
are absent and it is impossible to prove in such a simple way the convergence 
of iterations 

/0 rz 
y 2 exp((j>i(y))dyz = <p i+ i(z) + (z - x) 3 exp^x^dx 
-oo J — oo 

with 

00 = CX 

2.2 Uniqueness of solution for small arguments 

Now we shall prove the existence and uniqueness of non-iterated equation® 
(in the class of continuous functions) 

It can be done by several ways. One of the methods is to introduce the 
cut-off from the side of small z, i.e. to consider only z > —a with big and 
positive parameter a. It was done in pQ. Here we shall use more rigorous 
method. 

We shall show that for given c at 

1 c 4 

z < Zf in = -ln[— (1 - e)] 
with some small positive fixed e the solution of (0) exists and it is unique. 
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We construct iterations gi + i = P(gi)- The starting approximation is not 
important. We note that they do not go out of the class of positive continuous 
functions. Under the arbitrary initial approximation, all iterations starting 
from the first one will be positive. 

We have for approximations 

9i+i -9%= [ (x- zf exp(cx)(exp(-^) - exp(-^_ 1 ))rfx 

J — oo 

Then 

\9i+i-9%\< I (x - z) 3 exp(cx)\exp(-gi) - expi-g^ldx 

J -co 

It allows to write the inequality for norms ||...|| in the functional space 
C[— oo, oo] 

Hft+i - 9%\\ < I (x ~ z) 3 exp(cx)\\ exp(-^) - exp(-g- i _ 1 )||dx 

J — oo 

One can account that for positive gi and positive g i+1 the following in- 
equality 

I exp(^) - exp(g i+1 \ < \g t - g i+1 \ 
takes place. Having fulfilled the analogous transformations we come to 

||0*+i-0i||< / {x - z) 3 exp(cx)\\gi - g^i\\dx 

J —CO 

We take \ \g-i — <7i_i|| out of the integral and get 

\\gi+\-gi\\< I {x - zf exp(cx)dx\\gi - gi-i\\ 

J —CO 

Calculation of integral leads to 

6exp(cz) 

IWi+i -9i\\< 7 \\9i ~ 9i-i\\ 

c 4 

For the mentioned values of z we have 

\\9i+i-gi\\ < \\9i - e) 

The application of the recurrent procedure gives 

\\9i+i ~ 9i\ \ < (1 - e) 1 const 
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It is clear that the r.h.s. allows summation 

OO J 

y"](l — e) l const = -const 

Then this sequence is the Caushi sequence, then it converges, then the 
limit of iterations exists and it is the solution of equation (j3J) with given c. 
Since the iterations from the first one are positive, then the convergence takes 
place namely to this limit and the solution is the unique one. 

The proof of existence and the uniqueness of solution (jSJ) with given c in 
the region z < z in is completed. 

Now we shall consider solution at z > zi n . At first it will be difficult to 
consider all z, and we shall consider only those arguments, which correspond 
to the formation of essential part of the droplets. We shall determine them 
more accurately. 

Remark: We consider here a set of iterations. But this construction has 
no connection with a real rate of convergence of iterations with a fixed initial 
approximation. To estimate the real rate of convergence it is more correct 
to consider at first two different solutions, then to construct the iterations 
based on these solutions and to prove that they will come to one limit. Here 
the iterations are constructed correctly, but later the difficulty appears. The 
difficulty is that one can not use in constructions for calculations the limit 
of previous iterations as the base for the next type of iterations. One has to 
make a more detailed analysis including the establishing of the smoothness 
of the influence of the difference of the real iteration and the limit on the 
evolution at the further sub-period. It can be done but it is connected with 
some long formulas. So, we use the iterations, but keep in mind that we have 
to consider simply two different solutions. At the further sub-periods since 
the uniqueness at previous sub-period has been established one can choose 
the different solution with the common part at the previous sub-period. 

2.3 The properties of the function ip 

At first we shall present some properties of solution. One can notice that 

• The precise solution g increases if z increases 

• Any iteration gi increases if z increases 
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This follows from the explicit expressions. 

The mentioned properties allow to show some properties for the function 

ip = exp(c:r — g) and iterations ipi = exp(cx — gi(x)) 

• Function ip has only one maximum 

• Function ipi for every % has only one maximum 

• Function ip is always positive 

• Function ipi for every % is always positive 

• The maximum of ip is greater than the maximum of ip\ 

• The maximum of ip is less than the maximum of ip2 

Beside this one can show that if tp really exists, then all possible solutions 
lie between ipi and ip2- 

We shall prove that i/j 2 at big x has to decrease until zero. It can be seen 
from the following estimate. Certainly, 

/z 
(z — x) 3 ipi(x)dx) 
-oo 

for positive z is less than function 

ip2o = exp(cz — (z — x) 3 ipi(x)dx) 

J — oo 

The last function has the evident asymptote 

■^20 — exp(cz — z 3 / ipi(x)dx) 

which can be calculated exactly with some positive 

r° r° 6 
N + = / ipi(x)dx — / exp(cx exp(cx))dx 

J —oo J —oo C 

or 

JV+ = j[l-«p(-|)] 

So, we see that ip2 goes to zero rather fast and for given small positive 
fixed value of ip2 (let it be /) one can see the boundary Zf in where ip2 will be 
smaller than I. This completes the proof. 

Since one can prove that solutions ip of iterated equation lie below if) 2 the 
same conclusions take place for the iterated equation. 
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2.4 The uniqueness of solution for the intermediate 
values of argument 

Now we shall prove that for z in < z < Zf in the solution exists and it is the 
unique one. 

Here we construct the iterations of the special type 1 (marked by the 
subscript spl where it is important) as the following ones: 

• Before z < z in the initial approximation g is the precise solution (we 
already proved that it is unique). 

• At z > z in we take g = g(z in ) 

• The recurrent procedure remains the previous one. 

The old iterations will be marked by spO. 
These iterations satisfy the following properties 

• All iterations at z < Zi n coincide with precise g 

• Since for z > z in we have g > g(z in ), g sp2 > g(z in ) and the recur- 
rent procedure in the same one, we can come to the same chains of 
inequalities 



We denote by z maxi the maximum attained by ipi and by z max the maxi- 
mum of precise ip. One can see that 



g.s P i o < g.s P i 2 < ••• < g < ■■■ < g sp i 3 < g.s P i 1 



dg_ 

dz 




So, then 



max 




Z, 



■max spO 1 *^ Zi 



■max 



Z, 



■max 



< Z, 



and 



■max spl 1 ^ Z; 



■max spl 2 



Z, 



■max spl 1 *^ Z; 



■max 
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Zmax *^ Z max S pi 2 

For the difference gi + i — ^ one can write 

9i+i ~ 9i= I (z- xf exp(cx)[exp(-gi(x)) - exp(-^_i(x))]da; 

J— oo 

Then 

/z 
(z - x) 3 exp(cx){exp(-gi(x)) - exp(-^_i(x))]da; 
- in 

and 

/■z 

|ft+i-ft|< / (z - x) 3 exp(cx)\exp(-gi(x)) - exp(-^_ 1 (x))|dx 
Since > 0, > one can write 

|ft+i-ft|< / (* - x) 3 exp(cx)|5fj(a;) - gi_ x (x)\dx 

•> Zin 

and 

~9i\\ < / - x) 3 exp(cx)dx\\gi(x) - gi-i(x)\\ 

%in 

where the norm is taken in C[— 00, 00]. 

Since z — x < Zf in — z in one can come to 

/z 
exp(cx)dx\\gi(x) - gi-i(x)\\(z fin - z in f 
-in 

Since exp(cx) < exp(czf in ) we see that 

Ibi+i - 9i\\ < f dx\\gi(x) - gi_i{x)\\{zf in - z in f exp(cz fin ) 

and having calculated the integral we get 

\\9i+i-9i\ \ < (z - z in )\\gi{x) - gi-i(x)\\(z fin - z in ) 3 exp(cz fin ) 

The recurrent application of the last estimate (with the explicit integra- 
tion of (z — Zi n )) leads to 

II II ( ^ ^ITL ) / \ S / \ 

\\9i+l ~ 9i\\ < 7. COnst{Zf in - Zin) exp(czfi n ) 
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or simply to 



i\ 

Summation of the r.h.s. of the last relation leads to exp(z — Zi n ). Then 
this sequence in the Caushy sequence. Then it must have a limit. This limit 
will be the unique one and it is the unique precise solution of equation (jSJ). 

The uniqueness in the global sense can be proven by approach used in 
investigation of the iterated solution. It has to be repeated every time we 
come to analogous situation. 

The existence and the uniqueness of solution of equation © for z < zj in 
are proven. 

Generally speaking this proof is sufficient for uniqueness at every finite z 
but here we have estimated exp(cz) by exp(cz/j n . It is possible to give more 
precise estimate, which will be done below. 



\9i+i ~ 9i\ 



2.5 The uniqueness for the big values of argument 

Now we shall prove the existence and uniqueness for the rest z. 

We construct iterations of the special type 2. The procedure is absolutely 
analogous to the special type 1 but now the iterations are the precise solution 
until Zfi n . All properties mentioned for the special type 1 remain here. 

Now we rewrite equation for fa and have 

-fa+i + fa = ( (z- x) a [exp(fa) - exp(0j_i)]<ix 

J z fin 

Here a is the power 3. Now we keep parameter a because in further inves- 
tigations [2] it will be necessary to prove all for the arbitrary positive (and 
not integer) power. 

For the integer power the task is more simple because having differenti- 
ated several times we can kill the integral term and reduce the equation to the 
differential equation. Then we can use all standard theorems for differential 
equations. 

Then 

| - fa+i + fa\ < j (z - x) a \ exp(fa) - exp(0j_i)|ote 

•* Zfin 

Since at z > Zf in 4>i < fa < max fa = I one can write 

| exp(^) - exp(0;_i)| < r\fa - fa^\ 
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with some fixed r Then 

| - <f> i+1 + (f>i\ < ( (z - x) a \4>i - (f>i-i\rdx 

and for the norms in C 

|| - + < / (z - x) a dx\\(f)i - <pi-i\\r 

J Zfin 

Since (z — x) < (z — Zf in ) one can see that 



|| - + <j>i\\ < f dx(z - ZfinYWcpi - <pi-i\\r 

J Zfin 

Having calculated the integral we can finally come to 

|| - + <f)i\\<(z- z fin ) a+1 \\(t)i - 0j-i||r 

Having applied this estimate % times with explicit integration one can 
come to the following estimate 

( z _ Zf in ) w ^ 

\\<f>i+l -<f>i\\ < const 7T 

V{l) 

with some constant and two functions v and w. These functions have prop- 
erties 

v(i) > i\ 

w(i) < ([a] + l)i + consti 

where [a] is the minimal integer number greater than a and consti does not 
depend on i 

So, we can come to 

\\(f) i+1 -<j>i\ \ < const ^ 

Consider now z which satisfies condition 

Z - Z fin > 1 

Then 

( z _ z \{[a\+l)i+consti 

\\<pi + i - <pi\ \ < const 
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One can easily see that the term in the r.h.s. of the previous relation is 
the term in Taylor's serial for the function 

const(z - z fin ) M exp((z - z fm )^ +1 ) 

So, it is the Caushy sequence. So, the initial sequence is also the Caushy 
sequence. So, it converges and converges to the unique solution. So, the 
existence and the uniqueness of solution in the region z > Zf in + 1 is proven. 

One has also to note that in consideration of the second region there was 
absolutely no difference to prove the property until Zfi n or until Zfi n + 1. So, 
in the rest region Zf in < z < Zf in + 1 the existence and the uniqueness also 
take place. 



3 Existence of the root 

3.1 Some estimates for positions of the maximum of 
spectrum 

Now we shall give some estimates for position of maximum of ip at small and 
big values of c. 

Let z max be the point (for given c there will only one point) of the maxi- 
mum of precise ip. Let z max i be the maximum of fy. 
One can easily see that 

m(0) ° 



c 4 

Then the value z max \ can be found from the maximum of exp(ca; — ^ exp(crr)) 
and satisfies the equation 

6 ( \ 
c = —exp(cz maxl ) 

cr 

Then 4 

Zmaxi — In — 

c 6 

For every positive c it exists. 

The value of ipi at the maximum is 

Ipmaxl = 6Xp(cZ max i — exp( y CZ max i ) 
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and 



c 4 



In ip maxl = In — - 1 
6 



The last value goes to infinity when c goes to infinity. 
One can see that according to 

g'(z) = 3 / (z — x) 2 exp(cx — g(x))dx 

J —oo 

the following inequality 

g , (z)<g' 1 (z) 

takes place. One can easily calculate g[: 

rz Q 

g'i(z) = 3 / (z — x) 2 exp(cx)dx = — exp(c^) 

J-oo C 6 

From g'(z) < g[(z) it follows that 
At c = 6 we see that 

1 , 
Zmaxi = 7: In 6 > 
6 

So, at c = 6 we see that z max > 0. 

Since at every c the solution exists, since it is unique and since has one 
maximum then one can say that the dependence z max on c can be treated as 
a function. At c = 6 it is positive. 

3.2 The case of the small c 

Now we shall consider small positive c. We don't mention every time that c 
is positive but all conclusions for arbitrary c mean for arbitrary positive c. 
Since 

92 < 9 

one can see that 

Zmax2 ^ Zmax 

So, now we calculate z max2 . It is the root of equation 

/Zmax2 5 
(z ma x2 - x) exp(cx exp(cx))dx 
-00 <r 
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We introduce 

rz Q 

h(z) = 3 / (z — x) 2 exp(cx j exp(cx))dx 

J-oo ' c 4 

as a function of z. We see that I 2 is a growing function of z. The last 
equation can be rewritten as 

c = I 2 

If we change I 2 by some which is less than I 2 then the root z max2e of 
equation 

c = 7* 

will be greater than z max2 

Zmax2e ^ Zmax2 

and, thus, 

To construct we notice that 

exp(— — exp(cr)) > 1 — -^-exp(crr) 

since d 2 exp(x) / dx 2 > for all x. 
Then 

/•z 5 
I 2 — 3 / (z — a;) 2 exp(crr) exp( — - exp(cx))<ir 

J-oo C 4 

/■Z Q 

> 3 / (z — a;) 2 exp(cx)[l - exp(cx)]dx = I 2 

J-oo c 4 

Now we shall calculate the root of 1% — c. Having calculated integrals we 
come to 

6 6 6 

c = — exp(cr) — -^exp(2cx) — 

Having marked exp(cx) as y we come to the square algebraic equation 

4 a 36 2 

One can see that for small positive c there are two roots: one is very close 
to zero, another is greater than the first one but also it is rather close to zero. 
So, 
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Then 

exp(cx) ~ 

and cx goes to — oo. It means that z < z maX 2 e goes also to — oo. Since the 
roots to equation considered above are continuous functions of parameters 
(at positive c) we see that there exists some fixed small c y at which z maX 2 e is 
negative. Then at c y the value z max is also negative. 

If we now prove that the function z max (c) is a continuous function then 
according to Bolzano-Caushi theorem there will be a root z max = 0. 

4 The continuous character of the dependence 
of g on c 

4.1 The continuous character of the dependence of g 
on c at moderate arguments 

Now we shall prove that the solution of equation (J3J) depends on c continu- 
ously for c y < c < 6. Then it will follow that the maximum of cx — g will be 
continuous function and z max is continuous function of c. 

Imagine that g corresponds to c and g i corresponds to c'. Then 

g — g l = / (z — a;) 3 [exp(cx) exp(— g(x)) — exp(c') exp(— g l (x))]dx 

J — oo 

We see that the function 

q(c,g) = exp(cx - g(x)) 

at every x can be interpreted as a function of two arguments c and g. We 
see that the dependence on both arguments is the exponential one. Since the 
second derivative of exponent is always positive, one can give the estimate 
from above for q(c,g) — q(c c ,g c ). Namely, 

dq dq 
\q(c,g) -q(c',g c )\ < \{-Q-) m \\c- c l \ + \(-^-) m \\g - g'\ 

Here index m marks the maximal absolute values of derivatives which are 
attained on one of the ends of intervals [c, c'] and [g,g L ] correspondingly. 
Having calculated derivatives we come to 

\q(c,g)-q(c l ,g l )\ < exp(c*x) exp(-g*(x))\x\ |c-c'|+exp(c**x) exp(-g**(x))\g-g l 
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Here c* and c** are some values of c used in the maximal values of derivatives; 
g* and g** are some values of g used in the maximal values of derivatives. 
We don't require that g* = g(c*) and g** = g(c**). 
Then 

\d — 9 l \ < / (z — x) 3 [exp(c*x) exp(—g*(x))\x\\c — c l \ + 

J — oo 

exp(c**x) exp(— g**(x))\g — g L \]dx 
Having noticed that exp (—<?*) < 1 and exp (—<?**) < 1 one can come to 

\g — 9 l \ < f (z — x) 3 [exp(c*x)|x||c — c'l + exp(c**x)\g — g l \]dx 

Then for the norm in C we have 
1 1 C7 — ^7 ' 1 1 — / (z — x) 3 [exp(c*x)\x\\c — c l \ + exp(c**x)\\g — g l \\]dx 

J —oo 

or 

\\g-g'\\ < \c-c'\A + B\\g-g'\\ 

where z 

0<A= (z — x) 3 \x\ exp(c*x)dx < oo 

J — oo 

B= (z — x) 3 exp(c**x)dx 

J —oo 

The function B(z) is the growing function of z. We see that B = 
6 exp(c** z) I \c** 4 ) . The condition 

B = 1 - ei 

with a fixed small positive e\ determines the boundary Zj n *. then for z < Zi n * 
we see that ^ 

I \g — g l \ | < — \c — c'\const 
ei 

with some fixed const. So, g depends on c continuously at z < z ini( . 
One can also prove that 

x l exp(cx) exp(—g(x))dx i — 0,1,2,3 

-oo 

depend on c continuously. 

It can be done quite analogously. Certainly, the value of z in:¥ can be 
changed but remains the finite one. 
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4.2 The continuous character of the dependence of g 
on c at big arguments 

Now we shall investigate z > z in *. One can present g — g L in the following 
form 

\g-g'\=T 1 + T 2 
(z — x) I exp(cx — g(x)) — exp(c':r — g\x))\dx 

-oo 

T 2 — I (z — x) 3 \ exp(cx — g(x)) — exp(c l x — g l (x))\dx 

It is clear that T\ is the polynomial of the third order on z with coefficients 
proportional to U. Since coefficients are continuous functions of c one can 
state that the whole polynomial is the continuous function of c. So, one can 
write 

Xi < a±\c — c'| 

with some constant a\. 

Then for T 2 one can see that 

T 2 < f (z-x) 3 [exp(c*x)exp(-g*)\x\\c-c l \ + 

exp(c**x) exp(-g**(x))\g(x) - g l (x)\]dx 

Then since exp(— g*) < 1, exp(— g**) < 1 and exp(c*x) < exp(c*^), 
exp(c**x) < exp(c**z) one can come to 

T 2 < / (z — x) 3 [exp(c* z)\x\\c — c'| + exp(c**z)\g(x) — g l (x)\]dx 

Then since (z — x) 3 < (z — Zj n *) 3 one can come to 

T 2 < T 2 \ + T 22 

T 2 i — / (z — -2i„*) 3 |-2| exp(c*z)dx\c — c'| 

T 22 = j (z - z in *) 3 exp(c** z)\g(x) - g l (x)\dx 
Expression for T 21 and T 22 can be easy calculated analytically. 
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When we shall write the estimate for the norm in C. We simply have to 
write the norm \ \g — g l \\ instead of the absolute value \g — g l \ and then we 
can move \ \g — g l \ \ out of integral 

Then 

\\g-g l \\ <t 1 + t 21 + t 22 

where z 

T' 22 = f (z~z m ,) 3 ex P (c**z)dx\\g(x)-g i (x)\\ 

The integral in expression for T' 22 can be taken and 

T' 22 = (z-z m ^exp(c**z)\\g(x) - g l (x)\\/A 

Both Ti and T 2 i have the form 

const\c — c'| 

which is necessary and only T' 22 is expressed through \ \g(x) — g'(x)\\. So, to 
have the convergence after the recurrent use of the last formula one has to 
require 

(z - z in *f exp(c** x z) I '4 < 1 - e 2 

with a small fixed positive e 2 . 
This completes the procedure. 

We can fulfill this step again and again and z will grow. 
For any arbitrary fixed Zj iw we see that exp(c**^) < expfc**^^) and we 
have 

(^Z ^m*) ^ 



4 ( 1 ~^) u/4 



exp^c"2;/ in *; 

So, the size of the step at elementary procedure will be 

a, < ( ^.rV 

exp(c**Zf iw ) 

Since Az is finite, one can attain Zf in by the final number of steps. 

As the result the continuous character of dependence of g on c is proven. 
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4.3 The computational reasons for the uniqueness of 
the root of z max (c) 

Since the spectrum <fi has only one maximum, the continuous character of 
dependence of z max on c is proven. 

Since z max (c y ) < and z max (c = 6) > 0, the continuous function z max (c) 
must have a root. This root is a solution of the problem A. Thus, the existence 
of solution of the problem A is proven. 

But the uniqueness isn't proven. To do this we simply calculate z max as 
a function of c. This dependence is drawn in the Figure 1. 



1 

0.5 


/ c 

/ 2 4 6 8 10 







-0.5 




-1 
-1.5 




Figure 1: Coordinate z max as a function of c. 



One can see that the root is the unique one. The analytical proof can be 
found below but it requires another representation. 

One can see that for c < 3 the dependence z max (c) is a growing function. 
But for c > (Gexp(l)) 1 / 4 already the first iteration for spectrum has positive 
maximum. Since the maximum of precise solution can not be lower than the 
maximum of precise solution. The maximum of precise solution has to be 
positive. But the precise solution can not lie higher than ex. So, it must have 
maximum at positive x. It means that the slow decrease of z max (c) which is 
seen in the figure, can not lead to the negative z max . 
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So, the solution of the problem A is also the unique one. 



5 Iterated equation 

5.1 Some properties of solution of the iterated equa- 
tion 

For further purposes it will be necessary to investigate the iterated equation 
and to show the properties of solution. The iterated equation will be written 
in the following form 

/z rx 
(z - x) 3 exp(cx - / (x - yf exp(cy - g(y))dy)dx (9) 
-oo J — oo 

One can see the following property 

The solution of @) is so, that the spectrum, i.e. the function 

ijj = exp(cx — g) 

has only one maximum. 

Really, g is the integral with the positive sub-integral function. The 
function dg/dz can be presented as 



dg(z) 
dz 



z — x) 2 exp(cx — / (x — y) 3 exp(cy — g(y))dy)dx 



It is also the integral with the positive sub-integral function. Then the func- 
tion dg/dz is also positive. 
Now we recall functions 

2o = 

6 , , 



9i = ^exp^cxj 



92 



rz Q 

/ (z — x) 3 exp(cx exp(cx))dx 

J —oo C 



They are simply the already constructed iterations for eq. (J5|) . 
One can easily see that every solution of has to satisfy 

92 < 9 < 9i 

The proof is simply the comparison between expressions for g±, g2 and 
the iterated equation. 
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5.2 The uniqueness of the solution of the iterated equa- 
tion at small z 

Here one has to fulfill two procedures: 

f . To construct iterations to see the existence of solution 

2. To construct estimates for different solutions to see the uniqueness 

These procedures have similar technical realization. This is the reason why 
here we restrict the consideration only by the second case. The first case is 
identical in ideology to the procedure used for original (not iterated) equa- 
tion. There we solved only the first problem. The second case for the original 
(non iterated) solution can be solved as it is described here. 

Actually the existence for the iterated solution is already proven as the 
limit of increasing sequence (odd iterations) restricted from above and the 
limit of decreasing sequence (even iterations) restricted from below. 

The proof is the following: 

Suppose that g and q are two different solutions of equation ©. Then 



We shall formally denote the initial functions as go an d % an d determine 
the "formal iterations" gi and as 







cx — 



J —oo 




q i (y))dy)dx 



Then we have 





We consider z < z in . Recall that zt n is chosen from condition 



6 



j exp(cz in ) = 1 - e 



c 
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with a small positive e. 
We see that 

z in < z < x <y 
One can reorganize the previous equation as 

9i+\{z)-q i+1 {z) = [ (z-x) 3 exp(cx)[exp(- f (x-y) 3 exp(cy-cn(y))dy)- 

J — oo J — oo 

/X 
(x - y) 3 exp(cy - q { (y))dy)}dx 
-oo 

Then 

(^-a;) 3 exp(ca;)|exp(- / (x-y) 3 exp(cy-cji(y))dy)- 

-oo J —oo 

/X 
(x - y) 3 exp(cy - qi(y))dy)\dx 
-oo 

Since for positive a±, a<i one can see that 

| exp(— oti) — exp(— a 2 )| < \ot\ — CK2I 

and I-oo( x -y) 3e ^P( c y-9i(y))dy: /f^(x-y) 3 exp(q/-<fc(y))dy are evidently 
positive then 

(:r-y) 3 exp(cy-^(?/))cfr/)--exp(- / (x-y) 3 exp(c|/-g i (|/))cit/)| < 
-00 J— 00 

(x-y) 3 exp(cy-^(y))rfy) - / (x-y) 3 exp(q/-(Z;(y))dy)| 
-00 J —00 

Having rearranged the r.h.s. we come to 

/x rx 
(x-y) 3 exp(cy-gi(y))dy)-exp(- / (x-y) 3 exp(cy-qi(y))dy)\ < 
-00 J— 00 

(x - y) 3 exp(q/)| exp(-^(y)) - exp(-%(y))|dy 
Since qi(y) > 0, <&(y) > 0, one can see that 

I exp(-9i(v)) d v) - ex p(-Qi(y)) d y)\ < 
\Hy) -Qi(y)\ 
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and then 

/x rx 
(x-yf exp(cy-gi(y))dy)-exp(- / (x-yf exp(cy-qi(y))dy) \ < 
-oo J —oo 

/X 
(x - yf exp(cy)\gi(y) - q~i(y)\dy 
-oo 

Then for the norm in C space one can write 

/x rx 
(x-yf exp(cy-gi(y))dy)-exp(- / {x-yf exp(cy-qi(y))dy)\ < 
-oo J —oo 

/X 
(x - yf exp(cy)dy\\gi(y) -ft(y)|| 
-oo 

Having calculated the integral one comes to 

/X 
(x-yf exp(cy-qi(y))dy) \ < 
-oo 



I exp(- 

J —oo 

0- - e)\\Mv) - Qi(v)\\ 
Now we shall estimate — qi+i\. We have 

|ft+i-ft+i|< / {z - xf exp(cx)(l - e)\\gi(y) - qi(y)\\dx 

J —oo 

and then 

\9i+i - < (i - e) 2 ll^(y) - qi(y)\\ 

We can apply this recurrent procedure several times and come to 

\9i+i -qt+i\ < (i -t) 2l \\go(y) -qo(y)\\ 

So, we see that \gi+i — Qi+i\ can be estimated by some positive members of 
geometric progression. So, it means that this sequence converges and the 
limit is the real unique solution. The difference between ^ and & can be 
made negligibly small. 

So, we conclude that g and q is actually one solution. The uniqueness is 
proven. The existence can be proven by construction of Caushi sequence. 
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5.3 Further properties of solution 

We need some estimates for ip = exp(cx — g). 
Having constructed 

ii = CX- Cji 

i = cx-g 

we see that 

cx - gi < li < cx - g 2 
cx - gi < £ < cx - g 2 

Then it follows that 

max^i < max[cx — g 2 ) 

max^i > max[cx — g\) 

max I; < maxicx — g 2 ) 

max!; > max{cx — g\) 

and we have the estimates for the maximum. 
We introduce z 2 i as to satisfy 

(cx - g 2 ) | 221 = max{cx - gi} 

There are two z 2 \. We shall note them z 2 u and z 2 i r and choose to have 
Z211 < z 2 \ r . 

One can easily see that 

Z211 < Zmax i < Z 2 \ r 

where z max « is the point where the maximum of ^ is attained. It is also clear 
that 

z<mu < z max < z 2 i r 

where z max is the point where the maximum of £ is attained. 

So, now the natural boundary z 2 \ r appeared and we shall consider z < 

Z2lr- 

Since to solve the problem of existence we have to construct iterations 
then the estimates for iterations are presented. 



25 



5.4 Uniqueness for intermediate z 

Suppose that there are two different solutions g and q. According to the 
previous constructions they must coincide at z < z in . Then 

g(z) — q(z) — / (z — x) 3 [exp(cx — / (x — y) 3 



exp(cy - g(y))dy) - exp(cx - / (x - yf exp(cy - q(y))dy)]dx 

J —CO 

Then for the absolute values 

\g(z)-q(z)\< [ (z-x) 3 \exp(cx- [ (x - yf 

/X 
(x - y) 3 exp(cy - q(y))dy)\dx 
-co 

Then since z — x < z — z in we come to 

\g(z) - q(z)\ < (z - z in ) 3 \ exp(cx - / (x - y) 3 

J Z in J —CO 

/X 
(x - y) 3 exp(cy - q(y))dy)\dx 
-co 

The function 

/X PX 
(x-y) 3 exp(cy-g(y))dy)-exp(cx- / (x-y) 3 exp(cy-q(y))dy)\ 
-CO J —CO 

can be estimated as 

/X PX 
(x-y) 3 exp(cy-g(y))dy)-exp(cx- / (x-y) 3 exp(cy-q(y))dy)\ < 
-CO J —CO 

/X PX 
(x-y) 3 exp(cy-g(y))dy)-exp(- / (x-y) 3 exp(cy-q(y))dy)\ 
-co J —CO 

Since both 



and 



/X 
(x - y) 3 exp(cy - q(y))dy > 
-oo 

/X 
(x - y) 3 exp(cy - q(y))dy > 
-co 



are positive then 

/X PX 
(x-y) 3 exp(cy-g(y))dy)-exp(- / (x-y) 3 exp(cy-q(y))dy)\ < 
-co J —CO 
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/x rx 
(x - yf exp(cy - g{y))dy - (x - yf exp(cy - q{y))dy\ 
-co J—OO 



-co J — CO 

The last function allows the estimate 



/x rx 
(x - yf exp(q/ - g(y))dy - / (x - yf exp(cy - q{y))dy\ < 
-co J — CO 

/X 
(z 2 ir ~ z in f[exp(cy - g(y)) - exp(cy - q(y))]dy\ 
-co 

Since g and q are positive we come to 

/X 
(z2ir - z in f[exp(cy - g(y)) - exp(cy - q{y))]dy\ < 
-co 

/X 
[g(y) - q(y)]dy\ 
-co 

The final inequality will be 

/X 
[g(y) -q(y)}dy < ||g-^||(a;-^ n ) 
-co 

As the result we come to 

\g-q\< (z 2 ir - z in f exp(cz 2 ir)\\g -q\\(z- z in f f (x - z in )dx 

Having marked that (z — z in ) < (z 2 i r — z in ) and (x — z in ) < (z 2 i r — z in ) we 
see that 

\\g-q\ \ < (z 2 i r - z in f exp(cz 21r )\\g - q\\(z - z in ) 

Now we can formally organize sequential application of the last estimate 
and having applied the last relation i times with explicit integration of (z — 

^ we get 

\\g-q\\< (z 21r - z in f exp(cz 21r )\\g - q\\(z - z in ) l /i\ 
One can easily see that the last r.h.s. is the term in the Taylor's expansion 

of 

(z 21r - z in f exp(cz 21r )\\g - q\ \ exp(z - z in ) 

with all standard consequences mentioned above in such a situation. 
The uniqueness for z < z 2 i r is proven. 

In principle it is sufficient for all finite z 2 \ r . But below more fine estimates 
will be proven. 
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5.5 Uniqueness for the big values of arguments 

Here we shall use the function £ or if instead of cx — g. Again we mark by 
'"' the solution referred to the iterated equation. 

To prove uniqueness for z > z 2 \ r one can note that 

/X 
(x — y) 3 exp(^)dy) < exp(£ 2 ) < maxexp(^ 2 ) 
-oo 

Then it is clear that 

/x _ rx _ 

(x - yf exp(^)dy) < R\cx - (x - yf exp(£)dy\ 
-oo J —CO 

with some fixed constant R. This already solves all problems. But we shall 
present a detailed derivation. 

One can also note that the previous consideration will be valid not only 
for Z2\ r but for any finite z. For such a value z 21r we shall choose z or x at 
which already 

/X 
(x — y) 3 exp(^)dy < 
-oo 



-oo 

Since one can prove that cx — fl^x — y) 3 exp(^)dy goes to — oo at big z the 
last requirement could be satisfied. 

Let us suppose that there are two different solutions. For two different 
solutions £ and <p one can write 

\<p-£\< [ {z-x) 3 \{-cx+ f {x-y) 3 exp{fi)dy)- 

J —oo J —oo 

/X 
(x - y) 3 exp(£)dy)\dx 
-oo 

We can rewrite the last equation as 

_ rz rx _ 

\<p-£\< (z- x) 3 / (x- y) 3 \ exp(ip) - exp(£)\dydx 



Since at z < z 2 i r the solution is unique then ip = £ and we cut the lower 
limit of integration up to z 2 \ r - Then 

\<p-£\< [ {z-xf f {x -y) 3 \exp(fi) - exp(i)\dydx 
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Since it is known that both <p and £ have only one maximum and decrease 
at \z\ — > oo it is clear that both <p and £ are less than some constant and 
then 

I exp(^) - exp(£)| < Q\(p - || 
with some fixed 1 constant Q. Then 



and 

\\<f- 



-£\<Q f (z - x) 3 f (x - yf\$ - i\dydx 

J Z21r J %21r 

H\<Q [ Z (z-xf [ X (x-y) 3 dydx\\p-Z\\ 



The chain of inequalities can be prolonged 



\\<p-i\\<Q\ Z {z-xf f (z-y) 3 dydx\\<p-£\\ 

J Z'21r J z 21r 

The r.h.s. can be rewritten in a following way 

\\$ - i\\ < q( [' (z - x)*dx) 2 \\<p - i\\ 

Then since (z — x) < (z — z 2 i r ) one can see that 

\\<p-£\\<Q([ z (z - z 21r fdx) 2 \\<f - i\\ 

J Z2lr 

In other regimes of vapor consumption it is important to prove the unique- 
ness for the arbitrary positive power a instead of 3. So, we get the following 
inequality 

ii^ -in <q([ z (z-z 21r rdx) 2 \\^-£\\ 

J Z21r 

Having calculated the integral we come to 



\\^ - £\\ < Q(z - z 21r ) 2a+2 \\^ - £\\ 

Now we apply this inequality % times with explicit integration of (z — Z2i r ) 
and the result of this integration. Finally we get 

ll~ ell <>■ Dn( Z ~ Z 21r) W ^ I, ~ Z || 
\\<P-Z \\< RQ TTTTT ^0 -6 

V(i) 



1 The value of Q is close to 1. 
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where (fo and £o are initial approximations. For functions V and W one can 
get some estimates 

V(i) > i\ 

W(i) < (2[a] + 2)i + consh 

where [a] is the minimal integer number greater than a and consti is some 
constant which is independent on i. 

Now we can make the following remark: We can assume that earlier we 
proved the uniqueness until z 2 \ r + 1- Then z — z 2 \ r > 1. 

Then we come to the following estimate 

( z - z-, )(2[a]+2)i 

||^ - CII < constQ { - ^ {z - z 21r ) cM 

One can see that the r.h.s. is the member of the Taylor's serial for the 
function 

Qconst(z - z 21r ) constl exp((z - z 21r ) {2[a]+2) ) 

So, then the conclusion about the uniqueness and existence of solution can 
be also proven. 

6 /-representation 

One can rewrite the equation © by the substitution 

x — > cx 
z — > cz 

in a following form 

g(z) — f (z — a;) 3 exp(x — g(x))dx 

J —oo 

with / = c~ 4 . 

This representation is more convenient to show the uniqueness of depen- 
dence of z max on c. 

Suppose that the are two parameters /i and f 2 . Let it be f 2 > f\. 
Then 

rZ 

9fi( z ) = h {z- x) 3 exp(x - g h {x))dx 
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9h(z) = h I ( z ~ x ? ex P(^ - 9h(x))dx 

J — oo 

For gf 2 we construct iterations with initial approximation 

9h o(x) = 9h(x) 

and a recurrent procedure 

9h i+i( z ) = h I {z-xf exp(x - g h i(x))dx 

J —oo 

Then 

9h i(*) = hl ( z ~ x) 3 exp(x - g h (x))dx = (f 2 / ' h)g h {x) 

J —oo 

So, 

9h i(^) > 9h(x)(l + e) 

with finite positive small e. 

One can see the following important fact 
The following formula 

d \9h i -9f 2 ol 



> 



takes place. 

This is clear if one notes that 



/z 
(z - x) 3 exp(x - g h i(x))dx 
-oo 

grows when z grows. 
Really 

3(/2 - /i) / (2 - xf exp(x - g fl i(x))dx > 

J — oo 

This proves the inequality. 

One can also prove the following statement: 
There exists a point z\ n until which gj 2 2 > 9/2 
This follows from the explicit estimates for the first and the second iter- 
ations with zero initial approximations. It is known that at z — > —00 they 
come very close and estimate the solution very precise. Really, 

exp(c:r) 

9h -»■ 0/i as ~ 6/1 — - A — = 6/1 exp(x) 
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9h 9h as ~ 6/ 2 — ~ A — = 6/2 exp(x) 

Then 

9h 2 as = 9f 2 1 as = 9h as 

Since 

5/ 2 1 as > 9f 2 = 5/i 

we see that asymptotically 

Then the necessary point exists. This proves the statement. 
As the result we see that at extremely big negative z the following hier- 
archy takes place 

9f2 < <?/2 2 < ••• < 5/2 3 < 5/2 1 

So, all of them are between 5/20 and 5/21 at extremely big negative z. 
Later iterations can go away from this frames. We shall mark the coordinate 
when iteration number % attains the mentioned boundary by z\. 

One can easily see that either z\ goes to infinity or z\ exists. One can 
also see that 

z'i+l > z\ 

It is easy to see because the integration is going only for the values of 
argument smaller than the current one. This proves the estimate. 
A question appears whether the limit 

lim z\ 

exists. 

The existence of a limit means that at this limiting point (let it be z'oo) 
the limit of odd iterations differs from the limit of even iterations. But 
both the limit of odd iterations and the limit of even iterations belong to 
the solutions of the iterated equation. But we have already proved that the 
solution of the iterated equation is the unique one. So, the contradiction is 
evident. We come to the conclusion that there is no such a limit. 

The last conclusion particularly means that 

9h > 9h 
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Really, the analogous reasons show that the possible points of crossing 
9fi 2j+2, 9f! 2j+4 5 etc. with 2 j for every j cannot have a finite limit. So, 
for every finite z we see that 

9h > 9h 2j 

Since for every finite z there is such j which provides 

9h 2j > 9h o = 9h 

we see that 

9h > 9h 

This proves the necessary inequality. 

The last inequality leads to the important relation 



df 

for any given value of the argument. 



d9f >0 



7 The uniqueness of the root 

We analyze the question: how many solutions with different c (or /) can 
have the maximum of cx — g (or x — gf) at x = 0? 

The most evident answer is that there will be only one solution but this 
has to be proven. To prove it we shall use /-representation. 

The equation for the coordinate z max of the spectrum will be the following 

(z ma x - x) exp(x - g{x))dx 

-oo 

Then 

X rZmax 

/ = -(/ {z max -x) exp(x - g(x))dx) 

O J — oo 

Now we can calculate dj jdz max Then we shall take it at z max = and if 
we are able to prove that df /dz max < then the necessary property will be 
established. 

The possibility to take here z max = is ensured by existence of the 
root which was proven a few sections earlier. Now we see that the proof of 
existence of the root was really necessary. 
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For the derivative df /dz max one can get the following equation 



= _ o( / ( z max - x) exp(x - g(x))dx)~ Y 

dZ max " J —CO 

djl™ x (z max - xf exp{x - #(x))cfc 



c max <J J -co 

where 

y = 



max 



The direct calculation gives 

(z max -x)exp(x-g(x))dx- {z max -x) exp(x - g(x))—dx 

-co J —CO df 



In the last relation one has to consider as a function of x. 
One can also get for Y(0) the following equation 

F(0) = ~ l-co ^ 6XP( " - 

Now we shall calculate dg/df. One can get the following equation 

dcf f x f x dg 

— = (x - yf exp(y - g(y))dy - f (x - yf exp(y - g ( y ))—dy 

dj J -co J —CO dj 

One can rewrite the last relation as 

^ = g I 'f - f j_J,x - yf e*p(y - g(y))^dy 

Since dg/df > (it has been already proven) then the sub-integral function 
(x — y) 3 exp(y — g(y))jjjr is positive and the we have 

dg 



/ x tig 
(x- y f exp(y - g(y))-^dy > 



df 
Then 

Now we shall return to the function Y. One can rewrite it as 

/Zmax dg 
{zmax - x)[2 - {z max - x)—\ exp(x - g{x))dx 
-co dj 
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Then 

fO An 

Y(0) = (0 - x)[2 - (0 - exp(x - g{x))dx 

J-oo dj 

The sign of Y(0) is important. 
Then we consider the function 

We have to recall that x takes here only negative values. 
Then 

I > [2 + Xy] 

For p one can take the estimate 

9 < /exp(a;) 

going from the first iteration with initial zero approximation. Then 

I > [2 + a;exp(a;)] 

One can easily see that the function xexp(x) has only one minimum at 
x — — 1 where this function is — exp(— 1). 
Then 

I > 2 - exp(-l) > 

Then F(0) > and 

df/dz max \ Zmax=0 < 
It means that the uniqueness of the root is proven. 

8 Uniqueness of solution of the problem B 

Suppose that £ x and £ 2 are two different solutions of the problem 
B. Then (this is the statement) 

/ V 2 exp(-^(y))dy = / y 2 exp(-£ 2 (y))dy 

J —CO J —CO 

We shall prove this equality. 
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Suppose that these integrals do not equal, let it be 



•o 



/_ 



■0 



y 2 exp(-^(y))dy > 



y 2 exp(-&(y))dy 



— oo 



— oo 



Then for the problem A one can see that c\ > c 2 and fx < f 2 . Then (because 
there is only one root of equation z max (f) = 0) we come to the conclusion 
that at least one solution can not have maximum of x — g/(x) at z — 0. So, 
at least one solution is not a solution of the problem B. 

We come to the contradiction. 

This completes the proof. 

We can formulate another statement: 

The solution of the problem B is unique. 

It is clear because as it is proven 



has one fixed value. Then it is the problem of solution of equation (J3J) with 
fixed c. The last problem has the unique solution. 

This completes the justification of the uniqueness of solution of the prob- 
lem B. 

One can also see that the solution of the problem A is the solu- 
tion of the problem B. 

To see this one can simply differentiate the solution and then use the 
condition for derivative. 

Since the solution of the problem A is unique and solution of the problem 
B is unique one can easily see that: 

The solution of the problem A and the problem B coincide. 

9 Similarity of the forms of the size distribu- 
tions 

Having established the uniqueness of solution of equation (jSJ) with given c 
(or /), the uniqueness of the solution of the problem A and the uniqueness of 
the solution of the problem B we can investigate the correspondence between 
solutions with different c. 

The main result will be the following: 
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Actually all solutions of equation fl^D with different c (or /) are 
identical. 

Really, in equation (jSJ) we fulfill the substitution 

zc — > z 



and come to 




The next transformation will be the shift 

z -> z + A 
x — > x + A 

We choose 

A ~ 4 In c 

The factor (2 — x) can not change, the region of integration will be 
actually the same - from infinity up to a "current moment" z. 
Then the amplitude 1/c 4 can be easily canceled. 
So, the spectrum, i.e. the function 

$ = exp(x — g) 

has one and same form, which differs only by rescaling of the argument 
and the rescaling of the amplitude. The numerical simulations confirm this 
conclusion (see the figure 2) 

The shift in these pictures was introduced only in order to see three 
curves. When the shift will be zero then only one curve can be seen. The 
coincidence is practically ideal. 

The certain question whether there can be simply multiple repeating of 
solutions appears here. But the uniqueness of solution of equation (J3J) with 
given c and the property dg/df > established earlier ensure the absence of 
such repeating. 
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